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Preface 


The  purpose  of  this  study  was  to  determine  if  a specific  thermo- 
nuclear weapon,  detonated  at  such  a height  as  to  preclude  crater  forma- 
tion, would  leave  a signature  which  would  remain  distinguishable  for  up 
to  approximately  13  hours.  Primary  consideration  was  given  to  thermal 
ground  heating. 

I would  like  to  thank  Dr.  Charles  J.  Bridgman,  of  AFIT,  for  his 
support  and  guidance.  Major  Brian  G.  Stephan,  SAC/XPFS,  the  thesis 
sponsor,  is  due  thanks  for  his  support  and  invitation  to  discussions 
with  contractors  studying  the  DA/s  problem. 

I aa  also  indebted  to  Dr.  Gerhard  Rohringer,  of  the  General  Re- 
search Corporation,  who  supplied  me  with  material  which  supports  the 
accuracy  of  my  thermal  heating  model. 


Gordon  E.  Kelly 
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Abstract 

Prompt  thermal  heating  of  the  ground,  fission  product  deposition, 
and  neutron  activation  of  soil  were  studied  as  possible  indicators  of 
nuclear  weapon  damage  levels.  The  time  of  interest  was  the  period  of 
approximately  thirteen  hours  following  a 500  KT  low  altitude,  non- 
cratering burst. 

A computer  model  was  constructed  to  calculate  the  distribution  of 
weapon-induced  temperature  increases  within  the  earth.  Surface  thermal 
signatures  at  four  hours  were  detectable  and  ranged  up  to  Ai  K°  above 
ambient.  The  signature  had  decayed  to  background  levels  by  13  hours. 

Fission  product  heating  was  determined  to  be  unusable  as  an  indi- 
cator of  a weapon's  ground  zero.  Fission  products,  deposited  in  high 
concentrations  on  the  earth's  surface  as  from  severe  rain out,  may  alter 
the  early  thermal  signatures  through  beta  particle  induced  ground  heat- 
ing. Calculated  temperature  increases  for  severe  rainout  ranged  up  to 
35  K°. 

Neutron  induced  activity,  for  a given  soil,  was  calculated  using 

approximations  to  the  ORNL-4464  4rrr  scalar  flux.  Activation  of  a 

thin  layer  at  the  earth's  surface  was  assumed.  The  integrated  gamma 
24 

flux  from  Na  at  13  hrs  after  detonation  was  calculated  for  various 
altitudes  and  ground  distances  from  the  ground  zero.  Values  over 
ground  zero  ranged  from  8. 50  x 10^  Y's/cm2sec  at  an  altitude  of  100  ft 
to  1.02  x 104  at  1500  ft  altitude.  At  a ground  range  of  10,000  ft  the 
integrated  flux  was  32  Y's/cm  sec  at  an  altitude  of  100  ft.  A list  of 
fission  product  gammas  which  could  mask  or  confuse  this  signature  was 
compiled. 

*« 
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RESIDUAL  SIGNATURES 


PROM 

THERMONUCLEAR  AIR  BURSTS 
I.  Introduction 

Purpose 

This  study  was  performed  to  determine  if  a weapon  whose  height  of 
burst  (HOB)  precludes  crater  formation  leaves  a detectable  residual 
signature.  The  reference  weapon  is  a 500  KT  thermonuclear  device  hav- 
ing a 250  KT  fission  yield.  It  is  detonated  at  an  HOB  of  5000  ft  above 
ground  level.  Preference  was  given  to  the  possible  detection  of  ther- 
mal signatures  by  state-of-the-art  airborne  infrared  sensors.  However, 
the  study  was  not  restricted  from  suggesting  different  systems.  The 
time  frame  of  interest  for  residual  signatures  lasts  for  up  to  approx- 
imately 13  hrs. 

Definition  of  Terms 

In  this  report  "early  times"  refers  to  times  up  to  detonation  plus 
four  hours.  "Early-time  signatures"  refers  to  signatures  from  weapon- 
related  phenomena  at  approximately  four  hours.  "Late  times"  exist  from 
after  the  four-hour  point  until  detonation  plus  approximately  13  hrs. 
"Late-time  signatures"  are  those  at  approximately  the  13 -hr  point. 

"Pission  fragments"  are  isotopes  produced  by  fission  of  the  we&p- 
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on' 8 fissile  material.  The  term  "fission  products"  denotes  the  aggre- 
gation of  fission  fragments  and  all  subsequent  daughters. 

"Ground  zero"  is  that  point,  on  the  earth's  surface,  which  direct- 
ly underlies  the  burst  point. 

Background 

A damage  assessment/strike  (DA/s)  sortie  could  be  flown  against  a 
previously  attacked  target  as  insurance  of  target  destruction.  In  a 
general  nuclear  war,  a certain  number  of  aircraft  could  be  assigned 
missions  of  this  type  and  would  be  designated  DA/S  sorties. 

While  inbound  to  the  target,  the  crew  of  a DA/S  sortie  would  have 
to  determine  if  the  objective  had  been  destroyed.  If  the  crew  could 
determine  that  the  actual  ground  zero  (AGZ)  was  within  a certain  range 
of  the  target,  the  target  would  be  assumed  to  have  been  destroyed. 
Weapons  assigned  to  a target  lying  within  the  specified  radius  from  the 
AGZ  would  not  be  released.  If  the  target  was  outside  the  radius,  the 
target  would  be  struck  as  planned. 

In  some  cases,  targeting  constraints  might  require  that  certain 
allocated  weapons  be  detonated  at  an  HOB  wnich  does  not  allow  crater 
formation.  Such  conditions  might  not  produce  usable  visual  or  radar 
indications  of  the  weapon's  ground  zero.  If  DA/S  tactics  could  be  ap- 
plied to  these  cases,  destruction  could  be  assured  without  the  needless 
expenditure  of  valuable  weapons.  Acceptable  solutions  to  this  problem 
are  limited  by  the  practical  considerations  of  a DA/S  sortie  in  a gen- 
eral nuclear  war  environment.  A DA/S  aircraft  approaches  its  target  at 
high  speed  and  low  altitude.  The  decision  to  strike  or  withhold  weapon 
release  must  be  made  while  inbound  to  the  target  on  the  bomb  run. 
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Multiple  approaches  to  the  target  are  not  feasible.  Decision-making 
information  must  be  quickly  and  easily  obtainable  while  the  decision 
making  algorithm  is  kept  simple. 

Assumptions 

In  this  study  the  earth  was  assumed  to  be  flat  and  infinite  in 
depth  and  breadth.  It  is  thus  termed  semi -infinite.  It  was  assumed 
that  the  earth's  surface  absorbed  all  incident  pnenomena.  Absorption 
was  allowed  in  a thin  layer  at  the  earth's  surface.  Heat  was  trans- 
ferred into  the  earth  by  conduction  only.  The  earth's  thermo -physical 
properties  were  not  allowed  to  change  while  the  phenomena  under  study 
were  occurring. 

The  fireball  was  considered  to  be  a point  source  for  thermal  out- 
put calculations.  Firestorm  conflagrations  were  not  considered.  All 
heat  related  assumptions  were  made  so  as  to  maximize  the  absorption 
and  retention  of  heat  by  the  earth.  Prompt  thermal  heating  calcula- 
tions exclude  the  relatively  small  contributions  from  shock,  neutron, 
and  gamma  heating  of  the  surface. 
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II.  Prompt  Thermal  Heating  About  Ground  Zero 


A nuclear  weapon,  detonated  near  the  earth's  surface,  may  have  as 
much  as  of  its  total  energy  manifested  as  prompt  thermal  radiation 
(Ref  1»24).  This  radiation  arises  when  the  weapon's  soft  X-radiation 
(hv  S 20  keV)  is  absorbed  by  the  air  within  a few  feet  of  the  detona- 
tion point.  The  absorption  heats  the  surrounding  air  and  forms  the 
characteristic  fireball  which  radiates  at  several  thousand  °C  for  much 
of  its  life.  The  fireball  emits  characteristic  blackbody  radiation 
having  wavelengths  primarily  in  the  ultraviolet,  visible,  and  infrared 
portions  of  the  spectrum.  This  radiation,  emitted  by  the  fireball  in 
the  first  minute  or  less,  is  termed  prompt  thermal  radiation  (Ref  2s 2 6, 
75*31?) • The  total  prompt  thermal  radiation  or  thermal  fluence,  F, 
incident  on  a target  is 

F «=  VfT/%R2  (cal/cm2)  (l) 

where!  W = weapon's  explosive  yield  in  calories 

f = thermal  fraction  (fraction  of  the  yield  which 
manifests  Itself  as  prompt  thermal  radiation) 

T “ atmospheric  transmittance  (T  £ l) 

R «■  distance  from  burst  point  to  target  (cm). 

To  determine  if  a residual  signature  would  exist  at  ground  zero , the 
fluence  at  this  point  was  computed,  and  that  value  was  assumed  to  be 
uniformly  distributed  over  a semi -infinite  earth  such  that  F calories 
were  deposited  on  and  absorbed  by  each  square  cm  of  surface.  The  tem- 
perature distribution  in  time  and  space  may  be  calculated  from  a solu- 
tion to  the  general,  one-dimensional  heat  conduction  equation  (Ref  3tll)i 


4 


(2) 


arCx.t)  . k a^r(x,tl 

at  pc  dx2 

where!  T(x,t)  » temperature  in  °K  at  depth  x and  time  t 

k **  thermal  conductivity  (cal/cm  sec  °K) 
p = density  of  soil  (gm/ciP) 
c “ specific  heat  (cal/gm  °K). 

The  general  boundary  conditions  are 

- V-  «cv-  V -453%il]  . « 

L Jx  = 0 

where:  Q(0,t)  = net  rate  of  heat  flow  at  the  surface  (cal/cm~sec) 

\b  = rate  of  heat  absorption 
Qcv  " rate  of  convection  from  surface  to  air 
Q l “ rate  of  re-radiation  by  the  surface. 

At  great  depths  the  change  in  temperature  is  negligible.  Since  the 
deposition  of  heat  on  the  surface  is  assumed  to  be  uniform  and  infinite 
in  extent,  no  heat  flow  occurs  in  the  horizontal  plane. 

Analytic  Solution  to  the  Heat  Conduction  Equation 

As  a limiting  case,  the  total  fluence  was  assumes  to  have  been  de- 
livered as  a delta  function  in  time.  The  earth's  surface  was  taken  to 
be  insulated  with  the  exception  of  the  fluence  input.  Thus,  predicted 
temperatures  will  be  unrealistically  high.  With  these  limitations  Eq  (3) 
becomes 

Q<0, t)  - - kp^X.Oj  ^ . F6(t).  (4) 
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With  these  restrictions  Imposed,  Eq  (2)  Is  solvable  In  closed  form  by 
the  method  of  Laplace  transforms.  The  solution  Is 


T(x,t)  - £ /3St  .•x2/'rat  (5) 

where  a ° k/pc  Is  the  thermal  diffusivity  with  units  of  (cm2/sec) . By 
Eq  (l)  a 500  KT  weapon,  having  a thermal  fraction  of  0.43,  if  detonated 
at  5000  feet  above  ground  level  in  an  atmosphere  whose  transmittance  is 
0.95»  wiH  deposit  a total  of  700  cal/cm2  at  ground  zero.  Fig.  1 is  a 
plot  of  Eq  (5)  for  three  representative  types  of  soil.  Curve  I is  for 
sandy  clay,  which  is  soil  type  I.  Curve  V is  for  very  dry  soil,  and 
curve  VI  is  for  wet  mud,  which  are  soil  types  V and  VI  respectively. 

The  thermo-physical  properties  and  common  names  of  the  six  soil  types 
considered  in  this  study  are  listed  in  Appendix  A. 


Time  (sec) 

Fig.  1.  Insulated  Surface  Temperature 
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The  conduction  of  heat  into  the  earth  was  slow  enough  to  allow  high 
residual  surface  temperatures  to  exist  at  late  times. 

The  effects  of  convection  and  re-radiation  were  considered  next. 
The  earth's  surface  was  assumed  to  radiate  according  to  Stefan's  law 
(Ref  4:48): 

*rr  ■■  0€T^  (cal/cm2sec)  (6) 

where:  a *•  Stefan-Boltzman  constant  (cal/cm2sec  °K^) 

€ " emissivity  of  the  surface  (unitless) 

Ts  - surface  temperature  (°K) 

Free  convection,  assumed  to  take  place  at  the  surface,  was  taken  to  be 

Qcv  " hTs  (cal/cm2sec)  (?) 

where  h - 2.98  x iO-5  (cal/cm2sec  °K)  (Ref  5*180,444).  With  the 

addition  of  these  terms  to  the  surface  boundary  condition,  an  analytic 
solution  to  the  heat  conduction  equation  was  no  longer  possible.  The 
method  of  finite  differences  was  employed  to  obtain  solutions  at  a fi- 
nite number  of  points. 

Finite  Difference  Solution  to  the  Heat  Conduction  Equation 

A computer  solution  to  the  prompt  thermal  problem  also  allowed  the 
input  fluence  to  be  distributed  as  a flux  in  time.  The  standard  thermal 
fluence  spectrum  was  modified  to  allow  for  KX#  of  the  500  KT  weapon's 
prompt  thermal  yield  to  be  delivered  in  ten  seconds.  This  modified 
spectrum  and  the  flux  input  spectrum  are  contained  in  Appendix  B.  For 
the  case  under  study,  the  most  rapid  change  in  thermal  output  occurs  in 
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the  first  two  seconds  of  yield  generation.  The  largest  time  interval 
of  this  two  second  period,  over  which  the  change  in  output  was  nearly 
linear,  was  chosen  as  the  flux  input  time-interval.  It  was  found  that 
for  intervals  shorter  than  0.1  sec,  the  late-time  temperatures  did  not 
vary  significantly.  This  value  of  the  input  time  interval  required 
500,000  time  levels  to  extend  the  calculations  to  13  hrs,  A graded 
time  mesh,  having  a maximum  interval  of  five  seconds,  was  used  to 
shorten  the  calculation  time.  To  further  shorten  the  calculation  time 
a graded  spatial  mesh  was  employed.  It  consisted  of  35  nodes  with  Ax 
varying  from  0.2  cm,  near  the  surface,  to  10  cm  at  depths  beyond  50  cm. 
The  graded  space- time  grid  is  described  in  greater  detail  in  Appendix  C. 

The  heat  conduction  equation,  Eq  (2),  was  differenced  using  a sec- 
ond central  difference  approximation  to  the  spatial  derivative  and  a 
first  forward  difference  approximation  to  the  time  derivative.  Eq  (2) 
becomes  a four-point,  two-level  equation: 


Ti,**l  " Wi,  jTi+l,j+  (1  " Bli,j"  B2i,j)Ti,j+  ^l.jVl.J 


(8) 


where:  T.  . «■  temperature  of  the  (i, j)th  mesh  point 

l»  0 

01!, i " a0Atj/(AxiAx1^+  Ax^) 

^i,  j “ ***/ (AXI  + AVW 

- distance  from  the  (i-l)th  spatial  node  to  the  i^h 

Atj  - distance  between  time  node  j and  ( j + l). 

For  space-time  stability,  both  Sl^j  and  82^  j must  be  less  than  0.5. 

A complete  derivation  of  Eq  (8)  is  in  Appendix  D.  The  surface  boundary 
condition,  Eq  ( 3)1  in  differenced  form  is 
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Q(0,t)  - -k 


T*  j - Ti 


and  an  application  of  the  Newton-Raphson  technique,  described  in  Appen- 
dix D,  provides  the  surface  temperature  at  time  level  j: 


where t 


3C4T4+  C„(9.93  x 10-6)T1,33+  C„mx  + T„ 
1 D 2 — 1_Q 2 l 

4C.T3+  C (3.97  x 10‘5)T°*33+  1 
In  2 n 


Cj  - AXjCTCA 
C2  “ ^ 

Tn-1  “ (n+l)^  iteration  of  the  surface  temperature 
at  time  level  j 

TFLX  = weapon's  thermal  flux  input  into  the  earth  at 
time  level  j 

Tg  **  temperature  of  the  first  subsurface  node  at 
time  level  j 


The  computer  program  listed  in  Appendix  E was  assembled  to  calcu- 
late the  change  in  the  earth's  temperature  as  a function  of  time.  The 
program,  referred  to  as  HTE4,  can  be  used  to  calculate  temperature 
distributions  up  to  50,000  sec  (13.9  hrs)  after  detonation. 

The  following  results  are  the  temperature  increases  produced  in 
soil  types  I through  VI  by  the  500  KT  standard  weapon.  Table  I lists 
the  surface  temperatures,  rounded  to  the  nearest  half  degree,  at  four 
and  thirteen  hours  after  detonation  for  the  six  soil  types  listed  in 
Appendix  A. 
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Surface  Teaperature  (K°  above  ambient) 


fable  I 

Surface  Temperature  Increase  at  4 and  13  Hours 

Soil 

Temperature  (K°) 

Type 

4 hrs 

13  hrs 

I 

11.5 

5.0 

II 

7.0 

3.0 

III 

6.5 

2.5 

IV 

10.0 

4.5 

V 

4.5 

1.5 

VI 

10.0 

5.0 

Surface  temperature  increases  for  soil  types  I,  II,  and  V are 
displayed  in  Fig.  2 from  10  sec  to  1000  sec  after  detonation. 


Fig.  2 Surface  Temperature  Increase  from  10  sec  to  1000  sec 
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Figure  3 is  a plot  for  soil  types  I,  II,  and  V from  1000  sec  to 
50,000  sec.  Note  that  the  addition  of  convection  and  re-radiation  at 
the  surface  has  reduced  the  late-time  temperature  of  soil  type  V from 
the  highest  value  in  Fig.  1 to  the  lowest  in  Fig.  3. 


rig.  3.  Surface  Temperature  Increase  from  1000  to  50,000  sec 


Figure  4 illustrates  the  temperature  distribution  into  soil  type  I 
at  various  times  after  detonation.  The  three  plots  are  for  10  sec,  100 
sec,  and  1000  sec  elapsed  times.  The  effects  of  both  conduction  and 
surface  heat  removal  phenomena  are  evident. 
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Temperature  (K°  above  ambient) 


Depth  (cm) 


Fig#  4.  Distribution  of  Above  Ambient  Temperatures 
Into  Soil  Type  I 
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The  value  used  for  the  surface  emissivity  was  unity.  Actual  emis- 
sivity  values  may  range  from  0.8  to  0.95  depending  on  the  surface  cov- 
ering. The  use  of  a value  of  0.80  Instead  of  1.00,  in  HTE4,  increased 
both  early  and  late-time  temperatures  by  about  4#.  For  soil  type  VI 
this  was  the  equivalent  of  approximately  0.4  K°  at  four  hours  and 
0.2  K°  at  thirteen  hours.  For  the  same  soil,  a 20#  error  in  the 

I 

re-radiation  term  produced  temperature  differences  of  200  K°  during  the 
flux  input  period.  This  difference  decreased  to  less  than  0.2  K°  at 
late  times. 

Surface  temperature  differences  caused  by  convection  term  differ- 
ences of  20#  or  less  were  negligible  shortly  after  detonation  when 
re-radiation  dominates  the  heat  transfer  processes.  At  late  times,  a 
. 20#  difference  in  the  convection  term  produced  a calculated  surface 

temperature  which  differed  from  the  given  value  by  up  to  10#. 

An  error  in  the  surface  temperature  was  induced  by  the  calcula- 
tional  algorithm,  of  HTE4,  which  neglected  temperature  differences  in 
the  spatial  mesh  which  were  smaller  than  0.1  K°.  Decreasing  this  value 
to  0.01  K°  increased  the  surface  temperatures  at  50,000  sec  by  less  than 
1.0#.  Hie  following  inputs  resulted  in  late-time  temperature  differ- 
ences of  three  percent  or  lessi  a thermal  fraction  between  0.33  and 
0.43,  a value  of  the  atmospheric  transmittance  between  0.73  and  0.95, 
a burst  height  of  5000  + 500  ft,  and  an  input  fluence  error  of  + 20#. 

The  computer  model  was,  therefore,  relatively  unaffected  by  errors  of 
less  than  20#  in  input  parameters  or  surface  loss  terms. 

HTB4  was  used  to  model  the  Hiroshima  bomb  as  a check  on  the  code's 
accuracy.  The  weapon  had  a yield  of  approximately  20  KT  and  was  deto- 
nated at  an  altitude  of  18 50  ft.  The  flux  input  was  modeled  according 
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■ 


to  resolve. 

The  stability  of  this  model's  predictions)  coupled  with  Its  abil- 
ity to  predict  values  which  lie  within  those  found  at  Hiroshima,  indi- 
cate that  the  surface  temperatures  calculated  in  this  study  are  real- 
istic. The  calculated,  late-time,  above  ambient  temperature  increases 
were  generally  less  them  5 K°.  This  upper  limit  figure,  coupled  with 
the  calculated  background  temperature  variations,  indicated  that  ther- 
mally produced  signatures  would  have  decayed  to  background  level  by 
late  times.  The  existance  of  early- time  signatures  would  be  strongly 
dependent  on  soil  type  and  the  time  since  detonation. 


III.  Fission  Product  Heating 


By  definition,  no  surface-fireball  interaction  takes  place  for  an 
airburst.  However,  large  amounts  of  surface  debris  have  been  observed 
rising  into  and  being  ejected  from  the  top  of  airburst  fireballs.  This 
debris  does  not  become  highly  contaminated  by  fission  products  as  might 
be  expected.  Within  a few  seconds  after  detonation  the  fireball  takes 
on  a toroidal  shape  with  an  updraft  in  the  center  and  a downdraft  ar- 
round  the  outside  rim.  Surface  debris  which  is  entrained  upward  passes 
along  the  toroid's  central  axis.  Little  mixing  with  the  fission  prod- 
ucts contained  in  the  toroid  occurs.  Experimental  evidence  indicates 
that  surface  material  thus  entrained  contains  negligible  fallout 
(Ref  10 j 105). 

For  the  first  few  seconds  after  detonation  the  bulk  of  a weapon's 
fission  fragments  exist  as  a gas  or  vapor  due  to  the  high  core  temper- 
atures, T > 10®  °K.  These  gaseous  and  vaporized  atoms  are  carried  up 
by  the  rising  fireball.  As  the  fireball  cools,  the  fission  fragments 
nucleate  and  condense  into  particles  with  a maximum  diameter  of  about 
20  microns  (Ref  11 « 102). 

For  purposes  of  a limiting  calculation  (maximum  fallout),  it  was 
assumed  that  the  weapon's  fission  fragments  were  spherically  diverged 
from  the  burst  point,  and  a ground  zero  concentration  was  calculated. 
The  temperature  increase  due  to  beta  decay  wets  calculated  and  taken 
as  an  upper  limit.  It  was  decided  that  if  this  increase  was  negligi- 
ble, then  no  actual  signature  would  be  detectable. 
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Heat  Generation  by  Fission  Products 


The  rate  of  beta  decay  in  a gross  fission  product  mixture  follows 
the  Way-Vigner  1,2  law  (Ref  12*1 327) * 

■ a(t)  = AjT-1,2  (11) 

where  A^  is  the  activity  at  t = 1 and  r = (elapsed  time  at  A(r)/elapsed 
time  at  Aj).  Elapsed  times  are  measured  from  the  instant  of  produc- 
tion. This  relationship  is  a reasonable  approximation  for  times  great- 
er than  ten  seconds  after  fission.  In  this  development,  the  maximizing 
assumption  was  made  that  no  decay  takes  place  prior  to  10  sec.  Thus 
the  entire  production  of  fission  fragments  is  present  at  t = 10  sec 
(t  “ l)*  If  Eg  is  the  average  beta  energy  per  decay,  the  total  energy 
given  off  from  all  beta  decays  is 

*t«t  - v r v1,2  dT  - svi  <iz) 

where  time  is  measured  in  units  of  decaseconds. 

Ihe  total  delayed  beta  energy  from  the  fission  of  or  Pu^39 

is  approximately  0.044W*"  where  is  the  fission  yield  in  energy  units 
(Ref  13»ch  2,p  12).  The  activity  at  time  t is  then 

A(r)  - (0.044Wf/5^)r"1,2  (B’s/decasecond)  (13) 

Timofeev  and  Nesytov  (Ref  1^*-*  3)  give  the  8 activity  at  one  minute 
after  weapon  detonation  as 

o 

Ag  - 10  q (curies)  (l4) 

where  q is  the  TNT  equivalent  in  tons.  Values  for  this  estimate  and 
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those  calculated  by  Eq  (13)  are  within  an  order  of  magnitude  of  each 
other  for  values  of  between  0.3  and  1.0  MeV. 

The  weapon's  fission  fragments  were  assumed  to  be  deposited  on  the 
earth's  surface  prior  to  the  occurrence  of  any  disintegrations.  At 
T ■ 1 beta  decay  was  assumed  to  begin,  with  half  the  betas  being  ab- 
sorbed by  the  air  and  half  being  absorbed  in  a thin  layer  at  the 
earth's  surface.  Backscatter  from  both  media  was  neglected,  and  all 
beta  energy  was  assumed  to  manifest  itself  as  heat.  Beta  heating  of  the 
air  was  neglected. 

Spherical  Divergence  of  Fission  Fragments 

As  a limiting  case  of  the  ground-zero  signature,  the  weapon's 
fission  fragments  were  spherically  diverged  without  other  attenuation. 
The  ground-zero  concentration  was  used  to  calculate  the  increase  in 
temperature  due  to  fission-product  beta  absorption  by  a semi -infinite 
earth.  The  rate  at  which  0's  were  absorbed  by  the  surface  was 

Rg  B’  A(t)/8tt(H0B)^  (0's/cm2decasecond)  (15) 

The  thermal  flux  due  to  this  absorption  was  the  product  of  the  rate  of 
absorption  and  the  average  energy  per  beta: 

TFLX  « 3.50  x 10"4  Wfr‘1,2/(H0B)2  (cal/cm2sec)  (l6) 

The  computer  routine  used  to  calculate  the  temperature  increases 
may  be  found  in  Appendix  F.  The  earth's  temperature  increases  were 
calculated  in  a manner  similar  to  section  I.  Since  the  temperature 
gradients  encountered  in  time  and  space  were  not  severe,  a constant 
nodal-point  spacing  was  used  in  both  dimensions.  The  time  interval 
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used  was  five  sec,  and  the  spatial  interval  was  0,2  cm. 

The  differenced  heat  conduction  equation  for  a constant  space-time 
mesh  is 


iJ-H 


8rm.j+  “ 


»)Ti.j+eri-i.j 


wherei  8 " < l/2. 


(17) 


For  the  250  KT  fission  yield  weapon  detonated  at  5000  ft  HOB,  the 
peak  temperature  increase  in  soil  type  VI  was  37*3  K°  at  ten  seconds 
after  deposition.  Four  hrs  later  this  temperature  had  decreased  to  3.5 
K°  and  at  13  hrs  it  was  2.5  K°. 


Thermal  Signature  Alteration  by  Concentrated  Fallout  Depositions 

Early- time  thermal  signatures  may  be  altered  by  a weapon's  own 
fission  products  or  those  from  an  adjacent  target.  A weapon  detona- 
ted in  or  near  a severe  thunderstorm  could  have  its  fission  products 
subjected  to  high  percentages  of  rainout.  As  an  example,  consider 
that  all  fission  products  are  washed  out  over  an  area  equal  to  the  max- 
imum cross-sectional  area  of  the  fireball.  The  maximum  diameter  is 
(Ref  15*319) 

Dmax  " ^w°’4  (ft)  (18) 

where  V is  the  total  weapon  yield  in  KT,  Hie  resulting  thermal  flux 
from  beta  decay  at  the  earth's  surface  is 

TFLX  - 2.58  x 1<T12  v£r"1,2/W°*8  (cal/cm2sec)  (19) 
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I 


( 


where  W**  is  the  fission  yield  in  calories,  and  W is  the  total  yield 
c 

in  XT,  Calculations  were  made  using  the  computer  code  of  Appendix  F 
for  the  standard  weapon.  The  results  indicated  a temperature  increase 
in  soil  type  VI  of  35  K°  at  four  hrs.  For  a 50%  rainout,  the  increase 
was  19  K°  and  for  a 25%  rainout  it  was  10  K°.  These  figures  are  not 
appreciably  altered  by  the  deposition  time  so  long  as  it  occurs  with- 
in a few  minutes  of  detonation. 

Results  obtained  by  spherical  divergence  of  the  weapon's  fission 
fragments  were  negligible.  This  calculation  was  made  as  an  upper 
limit.  Therefore,  detectable  residual  signatures,  indicative  of  a 
weapon's  ground  zero,  would  not  exist  at  early  or  late  times.  Calcu- 
lations concerning  concentrated  fallout  or  rainout  depositions  indicate 
that  such  depositions  may  mask  or  alter  early-time,  thermally-produced , 
ground-zero  signatures. 
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IV,  Gamma  Activity  from  Neutron  Activated  Soil 


Neutron  fluence,  blast  overpressure,  and  other  prompt  weapon  ef- 
fects are  relatable  through  the  target's  range  from  the  detonation 
point.  The  activity  of  specific  soil  isotopes  formed  by  (n,y)  reac- 
tions is  relatable  to  neutron  fluence  if  the  soil  composition  is  known. 
Thus  a link  exists  between  induced  soil  activity  and  prompt  weapon  ef- 
fects. This  link  may  prove  useful  as  a DA/S  indicator. 

Fast  neutrons  are  rapidly  thermalized  in  soil,  with  peak  activity 
occurring  two  to  three  inches  below  the  surface.  The  isotope  selected 
as  an  indicator  of  activity  must  be  formed  in  an  abundance  which  per- 
mits measurement  at  DA/S  altitudes  at  late  times  after  formation.  The 
choice  must  be  made  based  on  photon  energy,  (n,y)  cross  section,  half- 
life,  and  concentration.  Lavrenchik  lists  the  seven  primary  soil  con- 
stituents shown  in  Table  II  (Ref  i6il8).  Some  elements  may  not  be 


Table  II 

Soil  Constituents 


Element 

Si 

Fe 

Ca 

Na 

K 

P 

Mn 

% by  Vt. 
(gm/cc) 

.74 

.13 

.10 

.075 

.07 

.003 

.002 

Atoms/cc 
of  Soil* 

i.522 

1.421 

1.521 

2.121 

l.l21 

6.019 

2.219 

Isotope 

Formed 

Si31 

Fe39 

Ca4* 

Na24 

K42 

p32 

Wb) 

0.4 

2.5 

0.6 

0.5 

1.2 

0.2 

13 

* 1.522  is  read  1.5  x 1022 


common  to  all  soil  types,  and  composition  will  vary  from  place  to  place 
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Ca4^  and  P^2  are  beta  emitters.  Silicon  comprises  25*75?  of  the 

30 

earth's  crust  by  weight,  but  only  three  percent  of  this  is  Si  which 
captures  a neutron  to  form  5i^  (Ref  17*B-135)*  This  nuclide  emits  a 
1.26  MeV  gamma  only  0.07%  of  the  time  and  has  a short  half-life  of  2.6 
hrs  (Ref  18:10).  Of  the  remaining  constituents,  Na2^  (n,y)  is  the 

best  indicator.  Sodium  is  only  the  sixth  most  abundant  element  world- 

23  24 

wide,  but  Na  is  the  naturally  occurring  isotope.  Na  has  a 15  hr 

half-life  and  emits  a 2.75  MeV  gamma  during  each  nuclear  transformation. 

If  only  those  reactions  indicated  in  Table  II  occur,  2.47%  of  the  inci- 

24 

dent  neutron  fluence  would  produce  Na  atoms. 

If  all  incident  neutrons  are  absorbed  in  a thin  layer  at  the 
earth's  surface,  the  activity  as  a function  of  distance  from  the  burst 
point,  R^,  and  time  would  be  (Ref  19*189) 

A(Rb,t)  *=  0.0247  F(Rb)\e"Xt  (Na24  y's/sec/cm2)  (20) 

where  X is  the  decay  constant  in  sec 

To  determine  the  neutron  fluence  at  a point  on  the  earth's  surface, 

a survival  fraction  was  calculated  from  0RNL-4464  (Ref  20: 316-317, 393- 

395)*  The  4rrr2  scalar  flux  was  calculated  at  1500,  1800  , 2400  , 3600, 

and  4800  meters.  From  the  values  of  the  flux  at  these  ranges,  the  fol- 

o 

lowing  approximation  to  the  4nr  scalar  flux  was  empirically  obtained: 

F = lOOe-4*^**0  ^b  (neutrons/source  neutron)  (2l) 

This  development  neglected  the  effect  of  the  presence  of  the  earth's 
surface  on  the  neutron  field  and  is  conservative  in  that  respect. 

Table  III  lists  the  values  obtained  from  ORNL-4464  and  those  calcula- 
ted from  Eq  (21).  It  can  be  seen  that  the  maximum  variation  is  14%. 
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Table  III 

O 

Comparison  of  4rrr  Scalar  Flux  and  F 

s 

in  (neutrons/source  neutron) 

Range 

Rb(m) 

4rrr2  scalar 
flux 

F 

s 

1500 

i.02"1 

9.35"2 

1800 

2.04"2 

2.32"2 

2400 

1.63"3 

1.42"3 

3600 

5.43"6 

5. 37"6 

The  scalar  flux  was  assumed  to  be  isotropic.  The  spherically  di- 
verged, surviving,  neutron  fluence  which  is  absorbed  by  the  ground  is 
(Ref  13«ch  2,p  50) 

9C  0^-4.65x10"^,  _ 

pn(Rv)  “ ** (neutrons/cm2)  (22) 

4ttR£ 

where  S is  the  total  number  of  neutrons  given  off  by  the  weapon. 

Flux  of  Na  Gammas  Encountered  by  a DA/S  Sortie 

The  2.75  MeV  gamma  flux  available  for  sensing  by  a DA/S  sortie 
would  be  the  spherically  diverged,  atmospherically  attenuated  activity 
integrated  over  all  space.  The  flux  at  the  point  denoting  the  aircraft 
location  is 

fy  " J /T[A(Rb,t)exp(-iiR J/Wh  d0R  dR  (23) 

' R^-e  QmQ  w • 

(Y's/c»2sec) 
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where i a " aircraft  altitude  in  cm  (ALT) 

- range  variable  of  integration  relative  to 
the  aircraft 

R^  = range  variable  of  integration  relative  to 
the  burst  point 

4 - total  attenuation  coefficient  for  the  2.75 
MeV  gamma  in  air. 

nay  be  written  in  terms  of  R&  and  other  quantities  as 

K ■ CRf“  alt2+  1)2 - 2(R2-  ALT2)1/^  cos(2tt  - e)]1/2  (24) 

where  D is  the  ground  distance  from  the  point  representing  aircraft  lo- 
cation to  the  ground  zero.  The  geometry  involved  in  this  problem  is 
diagrammed  in  Appendix  G.  This  function  was  Integrated  numerically 
using  the  AFIT  subroutine  SIMPD  which  uses  Simpson's  rule  to  approxi- 
mate a double  integral.  It  was  found  that  flux  values  were  not  sig- 
nificantly increased  by  increases  in  the  upper  limit  on  R beyond  a 
value  equal  to  (D2+  ALT2)^2+  HOB. 

Figure  5 is  a plot  of  the  gamma  flux  encountered  13  hrs  after  det- 
onation at  the  altitudes  indicatedi  for  ground  distances  up  to  10,000 
ft  from  the  ground  zero.  In  order  to  put  this  plot  in  perspective,  one 
must  realize  that  an  aircraft  traveling  at  a ground  speed  of  300  knots 
will  cover  the  10,000  ft  ground  distance  in  just  under  20  sec.  For  an 
aircraft  traveling  at  500  knots  the  transit  time  is  just  over  11  sec. 

In  a real-world  situation,  gamma-emitting  fallout  depositions 

1 

could  confuse  a sensor  if  the  fallout  pattern  were  similar  to  activa- 
tion patterns. 
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Extraneous  Gamma  Radiation  from  Fission  Products 

At  one  day  after  fission f the  total  gamma  activity  is  approximate- 
ly 6.6  x 10^  curies  per  kiloton  of  fission  yield  (Ref  l6* ** ***ll).  Radio- 

24 

nuclides  which  emit  photons  near  the  2.75  MeV  Na  gamma  were  examined. 

Table  IV  is  a compilation  of  the  fission  product  gammas  from  250  KT  of 

24 

fission.  Those  listed  have  energies  within  0.3  MeV  of  the  Na  gamma 
(Ref  18*32-82),  (Ref  21*106-108),  (Ref  22*87-92). 


Table  IV 

Potentially  Troublesome  Fission  Product  Gammas 

Nuclide  Half-life  y Energy  in  MeV  Rate  of  Occurrence  at 

and  (*  occurrence)  13  hrs  After  Fission 

(y's/sec) 


D 84 
Br 

fW». 

00 

• 

min 

if 

• 

CM 

(8%) 

3.70“ 

♦Kr87 

76.0 

min 

2.57 

(3  5*) 

2.0616 

-Rb88 

17.8 

min 

2.68 

(2.3# 

5.2816 

Rb89 

16.0 

min 

2.59 

(13*) 

2.575 

Y94 

20.0 

min 

2.57 

(1.5*3 

< 7.52* 

***Ag112 

3. 06 

(1.3*) 

< 6.583 

3.2 

hrs 

2.55 

(11*) 

4.0113 

Cs138 

32.2 

min 

2.63 

(9*) 

4.1212 

La140 

40.2 

hrs 

2.53 

(3*) 

2.1815 

La142 

92.0 

min 

2.55 

2.99 

s? 

7-80l6 

■>■55 

Pr1*6 

24.0 

min 

2.73 

(1.7*) 

2.099 

* Kr87  is  the  daughter  of  Br87  (56  sec) 

**  Rb88  is  the  daughter  of  Kr88  (2.8  hrs) 

***  Ag112  is  the  daughter  of  Pd112  (21  hrs) 

The  nuclides  listed  are  all  fission  fragments  except  for  the  three  in- 
dicated by  asterisks. 
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If  the  weapon  which  produced  these  activities  were  detonated  in 
or  near  severe  thunderstorms,  large  portions  of  the  fission  product 
inventories  would  be  deposited  on  the  earth's  surface.  For  purposes 
of  illustration,  a circular  area  of  deposition  was  taken  from  the  max- 
imum fireball  size,  Eq  (l8).  Rained-out  contaminants  were  assumed  to 
be  uniformly  distributed  over  this  area.  The  flux  encountered  by  an 
aircraft  at  an  altitude  of  ALT  would  be  (Ref  23 i 266) 

r2 

F « / [A  expC-pB  J/2RJ  dR  (25) 

v dad 

where i rl  - ALT 

r2  - (ALT^+  R2)1/2 
R ■»  radius  of  the  area  of  deposition 
A ” activity  per  unit  area 

- range  variable  of  integration  a3  in  Appendix  G 
4 - total  attenuation  coefficient  of  the  gamma  energy 
under  consideration. 

A value  of  5*81  x 10  ^ cm  * was  used  as  the  attenuation  coefficient  for 
all  gamma  energies  considered  (2.?5  + 0.3  MeV).  The  flux  as  a function 
of  altitude  was  integrated  numerically  with  the  constant  A separated. 
The  worst  case  was  taken  to  be  at  a point  over  the  center  of  the  cir- 
cular deposition.  The  results  are  listed  in  Table  V in  terms  of  the 
fractional  flux.  The  flux  at  any  altitude  is  the  product  of  the  frac- 
tional flux,  for  that  altitude,  and  the  surface  activity  per  unit  area. 

88 

Rb  is  potentially  the  most  troublesome  gamma  emitter.  For  a 
100%  fission  product  ralnout  from  a 500  KT  device  with  a 50%  fission 
yield,  the  activity  per  unit  area  Is  2.37  x 10^  photons/cm2sec  with 
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energy  of  2.68  MeV,  The  flux  as  seen  by  an  aircraft  over  ground  zero, 
from  Table  V,  would  range  from  1.57  x 10  to  2.35  x 10  gammas/cm  sec. 
This  rate  would  be  sufficiently  high  to  mask  signals  or  to  confuse 
logic  circuits  designed  to  act  on  rates  proportional  to  those  shown 
in  Fig.  5*  The  majority  of  such  gammas,  as  described  in  this  section, 
would  have  to  be  eliminated  by  selectively  measuring  the  2.75  MeV 
gamma  and  by  detector  characteristics. 


Table  V 

Fractional  Flux  of  Fission  Product  Gammas  at  Various  Altitudes 

Altitude  (ft) 

Fractional  Flux  (Fg/A) 

100 

.664 

200 

.393 

500 

.133 

1000 

.0333 

1500 

.0099 

Neutron  activation  of  soils  would  produce  detectable  late-time 
signatures.  Exact  knowledge  of  the  extent  of  soil  activation  or  the 
geometric  distribution  of  activity  would  be  impossible  to  obtain. 
Measurements  made  by  a gamma  detection  system  would  be  relative  to  a 
hypothetical  profile  as  in  Fig,  5*  The  sensing  system  would  have  to 
be  able  to  evaluate  all  possible  profiles  and  scan  the  vicinity  of 
the  target  area  to  determine  the  weapon's  ground  zero. 
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V,  Conclusions 


Signatures  from  Prompt  Thermal  Radiation 

The  prompt  thermal-heating  model  produced  calculated  temperature 
distributions  for  the  Hiroshima  bomb  which  are  consistent  with  histori- 
cally recorded  facts.  It  was,  therefore,  assumed  to  be  reasonably  ac- 
curate when  used  for  the  higher  yield  of  this  study. 

The  residual  ground-zero  temperature  increase  for  the  500  KT 
yield  was  a maximum  of  5 K°  at  13  hrs.  For  very  dry  soil  the  late- 
time temperature  increase  was  only  1.5  K°.  Since  the  heat  absorbed  by 
the  earth  could  reasonably  be  expected  to  have  a drying  effect,  late- 
time temperature  increases  would  be  in  the  vicinity  of  two  to  three 
degrees.  This  temperature  range  is  below  the  calculated  variation  in 
background  temperature  which  occurs  due  to  solar  heating. 

Since  the  assumptions  made  tended  to  maximize  heat  retention,  the 
values  obtained  must  be  considered  maxima.  Late-time  thermal  signa- 
tures would  not  be  adequate  or  reliable  as  DA/S  indicators.  Adequate 
signatures  would  exist  at  early  times,  but  the  possible  masking  effects 
of  thermal  heating  by  concentrated  fallout  depositions  from  nearby 
bursts  must  be  considered. 

Fission  Product  Heating 

The  assumption  of  spherical  divergence  used  in  the  calculation  of 
ground-zero  effects  was  unrealistic.  Most  fission  products  in  an  air- 
burst  remain  in  the  rising  fireball  and  subsequent  cloud.  In  spite  of 
this  assumption,  the  calculated  temperature  increases  about  ground  zero 
were  negligible  when  compared  with  the  calculated  background  variations. 


29 


Signatures  produced  by  concentrated  fission  fragment  fallout  or 
rainout  depositions  can  approximate  early-time  signatures  from  prompt 
thermal  heating  of  a weapon's  ground  zero.  This  would  be  especially 
true  if  the  deposition  were  graduated  to  heavier  concentrations  near 
the  center  of  a circular  area.  A location  where  surface  bursts  also 
occur  would  be  more  likely  to  experience  this  effect,  as  early-time 
fallout  from  airbursts  is  negligible. 

The  phenomenon  of  ground  heating  by  fission -product  fallout  leads 
one  to  consider  the  aerial  infra-red  reconnaissance  of  nuclear  battle- 
fields for  the  purpose  of  locating  heavy  fallout  depositions.  Since 
beta  and  gamma  activities  are  relatable,  the  heat  produced  by  beta  ac- 
tivity could  be  used  to  estimate  the  gamma  activity.  This  could  be 
used  to  estimate  the  potential  dose  to  personnel  entering  the  area. 

Gamma  Activity  from  Neutron  Activated  Soil 

Calculations  of  the  gamma  flux,  available  for  measurement  by  a 
DA/S  sortie,  indicate  that  usable  signatures  would  exist  at  late  times. 
These  signatures  would  be  measurable  at  altitudes  and  at  ranges  from 
which  DA/s  decisions  could  be  made  and  would  allow  cruise  missiles  to 
be  employed  in  a DA/S  role. 

The  system  to  detect  flu.  ] avels  would  have  to  be  gated  for  the 
24 

2.75  MeV  Na  gamma  to  avoid  the  extraneous  gammas  indicated  in  Table 
IV.  Due  to  the  short  reaction  time  available,  the  crew  indication 
should  be  a simple  go-nogo  indicator.  The  system  should  employ  a log- 
ic circuit  to  evaluate  cases  wherein  the  AGZ  lies  in  any  quadrant  rel- 
ative to  the  aircraft  position  at  decision  time. 

To  further  refine  these  calculations,  the  effect  of  the  earth  in 
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the  neutron  field  should  be  considered.  Gamma  flux  based  on  the  dis- 
tribution of  induced  activity  into  the  earth  should  be  calculated.  An 
exponential  atmosphere  and  an  attenuation  coefficient  which  varies  with 
it  should  be  employed.  The  addition  to  the  gamma  flux  by  scattered 
photons,  which  lose  only  negligible  portions  of  their  energy,  should  be 
considered. 
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Appendix  A 


Soil  Types 

The  soil  types  identified  in  this  thesis  were  chosen  to  offer  a 
wide  range  of  properties.  Some  types  may  occur  in  the  same  locale 
while  others  would  be  mutually  exclusive.  Table  VI  lists  the  six 
soil  types  with  their  thermo-physical  properties. 


Table  VI 

Thermo-physical  Properties  of  Soil  Types 

Soil 

Type 

Thermal 

Diffusivity 

a (cm  /sec) 

x 10"3 

Thermal 
Conductivity  k 

(cal/cm  sec  °K) 

x 10'3 

Specific 
Heat  c 

(cal/gm  °K) 

Density 

D 

(gm/cm3) 

I 

3.?o 

2.17 

0.33 

1.78 

II 

1.20 

.672 

pc  *= 

0.56 

III 

2.00 

.627 

0.19 

1.65 

IV 

1.90 

1.70 

0.53 

1.67 

V 

2.00 

0.40 

pc  *= 

0.20 

VI 

2.00 

2.21 

0.60 

1.50 

The  specific  heat  and  density  values  were  not  available  for  soil 
types  II  and  V.  They  are  listed  as  products.  Table  VII  lists  the 
common  names  and  the  references  from  which  the  six  soil  types  were 

taken. 
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Table  VII 

Soil  Type  Names  and  References 

Soil  Type 

Common  Name 

Reference 

I 

Sandy  Clay 
(15#  moisture) 

(24:189) 

II 

Clay land  Pasture 

(24:189) 

III 

Quartz  Sand 
(medium  fine,  dry) 

(24:189) 

IV 

Calcerous  Earth 
(4356  water) 

(6:288) 

V 

Soil 

(very  dry) 

(6:288) 

VI 

Wet  Mud 

(6:288) 

Appendix  B 

Thermal  Flux  Inputs 

The  standard  thermal  spectrum  relates  the  percent  of  total  thermal 

energy  emitted  to  the  scaled  elapsed  time  t/t  where  t is  the  time 

' max  max 

to  the  second  thermal  maximum.  This  value  for  a 500  KT  weapon,  detona- 
ted at  an  HOB  of  5000  ft  is  O.585  sec.  Fig.  6 shows  the  calculated 
curve  and  the  adjusted  curve  which  is  the  result  of  requiring  that  100# 
of  the  total  thermal  energy  be  emitted  in  ten  seconds,  as  was  the  case 
in  this  study. 


Time  Since  Detonation 


The  time  interval  between  flux  inputs  was  chosen  from  an  expanded 
version  of  Fig.  6 so  that  the  rate  of  change  from  one  input  to  another 
was  approximately  linear.  After  2.5  sec  of  yield  generation  the  weap- 
on *s  adjusted  output  curve  becomes  linear.  The  number  of  the  last 
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tiae  interval  before  this  linearization  was  denoted  LA.  Proa  tiae  in- 
terval LA  + 1 to  the  end  of  flux  input,  the  flux  Mas  decreased  at  each 
tiae  interval  by  a constant  amount  called  FINAGLE.  This  method  of  in- 
put was  used  to  decrease  the  number  of  inputs  to  the  computer  program 
in  Appendix  E.  Figure  7 is  the  resulting  vstrlation  of  input  flux  with 
tiae.  The  total  fluence  at  ground  zero  was  700  cal/cm  . 


0 2 4 6 8 10 


Elapsed  Time  (sec) 

Fig.  7.  Thermal  Flux  Input  Spectrum 
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Appendix  C 


Spatial  Mesh 

The  spatial  mesh  size  in  HTE4  was  varied  to  reduce  computational 
time.  Near  the  surface,  large  gradients  require  closer  spacing  for  ac- 
curate calculations.  At  greater  depths,  the  lower  temperature  gradi- 
ents allow  a coarser  mesh  without  loss  of  accuracy.  There  are  35  spa- 
tial nodes  per  tine  line  as  indicated  in  Table  VIII. 


Table  VIII 
Spatial  Meshing 

Ax 

* of 

Depth 

(cm) 

nodes 

(cm) 

0.2 

11 

0-2 

0.5 

6 

2-5 

1.0 

5 

5 - 10 

5.0 

8 

10-50 

10.0 

5 

50-100 

Time  Mesh 

The  time  spacing  for  the  first  ten  seconds  was  kept  variable  to 
allow  the  flux  change  from  one  tine  level  to  the  next  to  be  kept  ap- 
proximately linear.  The  first  time  level  was  taken  to  be  t ■ 0 so  that 
for  Atj  “ 0.1  sec  101  nodes  were  required  in  the  first  ten  seconds.  If 
Atj  ■ 0.05  sec,  201  nodes  sure  required.  Table  IX  lists  the  tine  node 
spacing. 
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Table  IX 
Tine  Meshing 


Tine 

Increment 

Indicator 

Tine 

Increment 

(sec) 

Number 
of  Nodes 

Time  Period 
(sec) 

Atl 

variable 

(lO/dt^-H 

0-10 

< 

1.0 

90 

10  - 100 

“3 

2.0 

450 

100  - 1000 

5.0 

9800 

1000  - 50000 

Appendix  D 


Derivation  of  the  Differenced  Heat  Conduction  Equation  and  the  Newton- 
Raphson  Treatment  of  the  Surface  Temperature 

The  central  difference  operator,  6,  is  defined  by  (Ref  25*56) 

6y(x)  - y(x  + h/2)  - y(x  - h/2)  (26) 


Let»  x «*  iAx,  t - jAt,  then  x + h/2  ■ Ax(i  + i/2).  In  this  notation 
(time  dependence  suppressed) 

6Ti  " Tm/2  " Ti-l/2  to) 

If  Ax  is  not  constant,  the  central  difference  approximation  to  the  first 
spatial  derivative  is 


»T  „ fin/2  ' Ti-l/2 

dx  (Ax±+  Axi+1)/2 


where  the  Ax^  are* 


AXi  Axi41 

— v * *. 

i-i  i m 


(28) 


The  central  difference  approximation  to  the  second  derivative  is 


^~2  ~ Ax±+  Ax1^1^6Tl-^/2^Axi+l^"^6Ti-i/2^Axi^  ^ 


ix 


Carrying  out  the  indicated  operations  and  resurrecting  the  suppressed 
time  dependence,  the  above  equation  becomes 


i^t)  2(Ti+i.r  ru.\)  2(Ti,r  Vi.j* 

Ahd*Ui+  *4«  *4+  AV*!-* 


(30) 
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In  general,  a difference  equation  must  not  be  of  an  order  higher 
than  the  derivative  It  replaces  If  Instabilities  are  to  be  avoided. 
Thus  a first  forward  difference  operator  in  time  was  chosen,  and  the 
differenced  form  of  the  heat  conduction  equation,  Eq  (2)  is 


.J*r  T1..1 

At 


(31) 


where  At^  is  the  time  interval  between  nodes  j and  j+1,  and  a is  defined 
to  be  k/pc.  The  temperature,  for  points  below  the  surface,  at  time  line 
j+1  in  terms  of  the  temperature  at  time  line  j is  a four-point,  two- 
level  equation: 


2aAt, 

where:  81.  . J 5—  < l/2 

J Axi+1 


82, 


2oAt 


1 


i,J  A*i+<Vxi+i 


< 1/2. 


Hie  difference  approximation  to  the  surface  boundary  condition, 
Bq  (3)  is 


- \b-  v ^ 


(33) 


where  is  the  thermal  flux  absorbed  at  time  level  j.  Taking  all 

temperatures  to  be  relative  to  T.  . and  using  Bqs  (3),  (6),  and  (7), 

i»l 
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the  equation  for  the  surface  temperature  Is 


Ti,j  ■ 2-?8  * <»'5  *i;f) + t2,j  <») 

Since  the  surface  temperatures  encountered  are  large , the  surface  tem- 
perature at  time  level  j cannot  be  used  to  calculate  the  reradiated  flux 
at  time  level  j+1 . The  result  of  this  method  is  unstable  and  leads  to 
negative  temperatures.  The  Newton-Raphson  technique  was  used  to  avoid 
this  instability  (Ref  26:308). 

In  this  method,  a sliding  tangent  is  used  to  find  the  zero  of  a 
function.  The  following  equations  of  the  function  f refer  to  Fig.  8. 


(35) 

(36) 


Fig.  8.  Sliding  Tangent 


:’H>  ■ 


*•(0^)  - o 

V %n 


then 
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Appendix  E 


Computer  Program  HTE4 

The  basic  calculation  of  temperature  in  HTE4  involves  the  follow- 
ing! 

1.  Calculating  the  values  of  Rl.  . and  R2.  . for  the 

1 • J 1 » j 

time  increment  in  use 

2.  Calculating  the  subsurface  temperatures  using  the 
differenced  heat  conduction  equation 

3.  Calculation  of  the  surface  temperature  based  on 
flux  input,  convective  losses,  radiative  losses, 
and  the  value  of  the  temperature,  computed  in  step 
2,  of  the  first  subsurface  node 

4.  Iteration  of  the  surface  temperature  to  convergence. 

Throughout  the  program,  J is  the  time  index  and  I is  the  spatial 
index.  The  parameters  C^  and  C^,  defined  in  Appendix  D,  were  used  to 
reduce  the  physical  space  required  to  write  out  the  surface  tempera- 
ture equation.  Flux  input  begins  at  time  level  two,  and  the  choice  of 
1000  °K  as  a beginning  temperature  for  iteration  is  arbitrary. 

The  program  was  written  in  Fortran  for  use  on  a CDC  6600  com- 
puter. 
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SET 

TNE7  = TJ(l) 
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— comment — 
the  following  section 
calculates  temperature 
during  the  remainder  of 
the  flux  input 


BETA2(l  i 
GE.|5 


/ PRINT 
/BETA2(I)  GE..5/ 
/for  DELTAT(2)  / 


^ TEST  x 
JETA1/2(I] 


LT.  .5 


bstai(i) 
GE..5  1 


/ PRINT 

/betai(i)  ce.  .5 
for  DELTAT(2) 


J=NTFLX1 

N1 


(b  (b 


1 

r 

DO 

1 - 

ND. 

4 

190 

2. 

X 

CALCULATE 

TJl(l) 

subsurface 

temperatures 


temperature  calculations  from  100 
to  1000  and  from  1000  to  50,000  sec 
are  similar  to  the  proceeding  sec- 
tion. DELTAT(3),  DKLTAT(4) , MDTC, 
and  NDTD  are  used  to  calculate  BETA. 
The  remainder  of  the  program  is  out- 
put control  for  the  temperature  read- 
outs. As  it  is  listed,  it  is  valid 
for  DELTAT(l)  of  0.1  or  0.05  sec. 


3L0?*A?Y  Op  USED  IN  T»J_  PROGRA M MTE4-CYCLE  *5 


ALpMA  - TH"omAL  nirr||STVlTv  (CM**?/src> 

HiTSl(T)  _ TT  MENEIONLE??  PA^AM^Tr^s  TM  THE  MEAT 
B-T9?(T)  _ RONDIPTTON  FOUATTOM,  THEY  APE  FUNCTIONS 
0^  OELTAX(I),  OELr  AY('M>  , ALPHA, 

ANO  DELTAT(J). 

DELTA  T 1 - IT*"  TNT^PVAl  (S~n  OU^INS  FLUX  INPUT 
DELTA  T? 

DELTATY  - rr*r  IfJC °E MINTS  (SrC)  0p  IMOOEASINO  ST7F 
DELTAT4  JHIOH  A=»r  USP  0 TO  ^CREASE  COMPUTATION  TIMf 

rp5TLOM  - THE  EMTSSTVITY  OF  'HE  SURF ACE 

FINAGLE  - LINEAR  FLUX  r MANGE  IN~9FMEHT 

DELTAX(I)  - -‘▼STANCE;  Or  THr  r-TH  SPATIAL  NODF  FROM 
THE  OProrEOTNG  Mon? 

K - THERMAL  0 ONOUCtI VIt Y (CA. /EM  SEO-OFO  <) 

LA  - LAST  TH'OHAL  rLUX  INPUT  n^FCRE  INPUT  eLUX  BECOME? 
LINEAR 

N1-M1* 

Ml-M?  - COUNTER*?  USED  TO  ppI'IT  SURFACE  TTMpeoATURES 
AT  ▼-•-  INOTCATEO  POINT?  IN  TI*£ 

NDTA  - MUMRPR  OT  r IMF  TNTEPVAL?  OF  SIZE  DELTA T1 

NTT?  - NJMBE’  OT  TINE  IMTE®V1LS  OF  SI7E  DELTAT? 

NOTC  - NUN’pp  OT  TTME  TNTrRVALS  DF  ST*F  OELTAT3 

<nTI)  - NUMDF * OT  TTMF  TNTFpV«Ls  *)p  ST’E  DELTAT4 

NDX  - NUMD"R  Op  OFOTM  INTERVAL' 
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I 


x - N.i<vr  op  *mpr»«ai  pljy  impute 

S*  - PROOTT  Op  THC  OTppAM  '»HT7*44M  PONET4 NT  ANO  THr 

pmic!;i  vt”y  op  TM-  pn-opa'p  (C4L/c,'**R-o::o-''cr.  <**«»» 

T0C  - TIFrror^P  n-H Mr  EN  SMO'PSSIVE  TTE’MTTONS  OP  THr 
SU=>P40:  TrMP"R«TIJO~,  ',rLOW  W«ICM  ITERATIONS 
49 r f f?htM4 TPO 

TPLX(J)  - TH"°MAL  rL'JY  4T  TTH?  LIN;  J <CAL/0M**2-S"n 

T J ( T ) - AN  APP4Y  PO3  tHe  TEM3P0ft  RY  PTORAOE  Oe 
TM-  *ALUr«  OP  T(I,J>  ft T TIM:  LTN“  J 

TJ1(T)  - 4M  4P?ft Y P0»  TM"  Tr^nrjoR^Y  PT0R4SE  OP 

TM"  VALUES  OP  T(i,.JH)  rt  TIME  LEVEL  J*1 

TMIN  - S'JpPftOP  t"MPpP4TU^p  Or%<)  OPLOM  MHIOH  NO  PJ3TM£R 
CALOUl  4~I0NS  ARP  Mftnr  c 1 , T 3 0 - *»3, 1 10  SEC) 

TO  - «?'|PP4  3P  T-MPPoftTiJ3"  AT  TT*E  ZE°0  ( "OUAL  TO  C OES-O 

TS(J)  - TM£  PtP.PACP  TE^PPRA T’JRT  4T  TIME  LEVEL  J 
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Program  Listing  HTE4  (Grouped  according  to  the  operation  performed) 


R^A:  < 

OT^NSTON  TJ<35) , T J1 < 15) , *~TA 1 < ?5) , 0ETA2 < 35) , 
GOrLTAK(35)»T<>(n550)  ,TPLX  (?•“!) 

RFA0%  AlOHAtK,£OSU0M 

RrA0*,f3:TLTATl»,)'*LTAT2»DELTATl,3ELTAT4,N0TA,  NOTB , NOTC » NOT'D 

RFAO*,MOX,TOe,T'1IM 

RFf‘'*.MTFLX,LA*FI,JA'»Lr 

R^"**,  (TPLX  (J)  , J=lf L A) 

NP*1  r'JPX^l 

Rr«o*,  (PTLTAXd)  ,I  = l,N0Xi) 

LA1  = LA4-1 

00  I C J=  L Ai * MTpLX 
TPL< (J) =TFLX( J-i) -FINAGLE 
1?  GOMrTMJ* 

TTl=NTrA*0"LTATl 
TT?r NOTO*OELT*T? 

TTls  MnTC*  nELT  ATI 
TTUtH3T0*0: LTAT4 
TTJ?  eTTl  + TT? 

TTl?3=TTlfTT?*rT3 

TT=rTl*TT2fTT3tTT4 

P*I'IT* 

PPDT*,  "ALPHA  = “.ALPHA,” 

PPT'IT»,”K  = ”,K,”  C*L/$r C*CM*OES  <" 

PPIN  T * , “ NMMR*  9 OF  Or  PrH  INCREMENTS  = ",NOX 
PPI'fT*,“TIHf  INCREMENTS*” 

PPI'iT*,  "DELTA!  l 2 " f DEL  TAT 1,  “ SEC.  FROM  0 TO  ". 
CTTl,“3E:.“ 

PRINT*, "OFLTAT  2 = ", DELTAS, “ SEC.  FROM  ",TT1, 

C”  TO  ",TT12,”  SEC.” 

PRI')T»,"0£LTAT  3 = ”,0ELTAT3,”  SEC.  FPOM  ”,TT12, 

C“  TO  “,TT1’3,“  SEC,” 

poIv|T»,"r)^L.TAT  4 = ”, DELTA  T4,  ” SEC.  FROM  ”,TT123, 

C”  TO  ",TT,"  SEC.” 

LA3=LA*3 

PRINT*, "THE  FJORT  ",LA3t"  FL'JX  IMPJTS” 

PPT'fT  15,  (TFLX(J)  .J-1.LA3) 

15  FORMATES. 2) 

PRINT*, “SPATIAL  MESHt” 

«»PTN*»,  (OELTAX(I)  ,Isl,NOXl) 

PPINT* 

SE=(1.35E-i?) *FPSILON 
T0=?.  . . 

00  T *i »M0X1 
TJ(T )=T0 
TJ1(T)*T0 
20  CONTINUE 

TS (t  )*T0 

C-— — TrM»r»ATl)RE  CALCULATION  AT  SECOND  TIME  LEVEL—  — — 
C1»0FLTAX (?)*SE/K 
C?«0ELTAX(2)/K 
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TJ(i)  = (3*Cl*(Kni.**4)  *C2*(°.  93E-J6)  * (1 0 3 0 . *•  1. 33  3) 
C*C2*TFLX(2) 

C*TJ(2))  /(4*riM13,'C.**3)*C2*(3.97E-0  5)*(10  0 0,**.  3 33>*l.) 
?5  TH- >♦  = T J C i > 

TJC  )=  C**Cl*  (TNEW  **4)*C2*(''.<nE-06)*(TNEW  ♦*1.  333) 
C*C2,TFL<  (?) 

C*T J( 2) > / (“*Cl* (TNEW  **?)  + C2*(  3.  97E-15)*  (TNEW  **.  333)*1) 
IF(in^(TJ(l)-TM£K).SE.TOP)  GO  To  25 
T$<>>  = TJ(1) 

C Tr“3r  RA  T»JRr  CALCULATION  CURING  FLUX  INPUT 

RETAKl)  = C. 

«*FT!  ?(1)  =0. 

no  7c  i=?,mox 

R-m(I)  =?*AL3HA*OELTATl/(DCLTAX  (I)  *DELT4X  (1*1)  * 

COrLT  AX { I * 1 ) ** ?) 

IF  nr-TAi  (I)  ,r,-.,5)  GO  TO  29C 

ntT12(I)  =?»  ALPH\*‘r>ELTATl7((nELTAX(I>  **2)  fOELTAX (I) 
c*nr.TAx(i*D) 

IF(’,ETA2(I)  • GF. ,5 ) no  T0  291 
3D  CONTINUE 

00  *0  J=3,NTfLX 
00  ?T  1=2, NO* 

TJKDsPETAl  (I)*T  J(I*l)*(l-«»-TAl  (I) -9FT A2( I) > *TU( I) 

C*nrr6?(r)  «TJ(T-n 

IF  (AR5(TJ1(T) -TJ1 (I*l))«LT..i)  GO  TO  60 
50  CONTINUE 

TJK  N0X1)  =TJ1  (N)X) 

63  TJ1( 1) = ( 3*C1* (TJ( 1)  * * 4) *C2* <9. 9 3E-0 6 ) * (T J ( 1 ) ** 1 . 3 33) 

C+02*  TFL  X ( J) 

C+TJ1  (2)  ) /(4*Cl*  (T  J (1 ) **3)  *22*  (3.  O^E- 05)  * ( T J (1 ) **.  333)  *1) 
65  TNF4sTjkii 

TJK  1)=(3*C1*  (TN*W  **4>*C2*  (9 , 9 3E-0 6) * ( TNEW  **1.333) 

C*r  2*  TrL  X ( J) 

C*T.J;  (?) ) / (4*01* ( T NFW  **3)*C?*(3.97E-05)*(TNEW  **.  333)  *1) 
IF  ( A9?(TJ1 (1)-TNFW) .GE.TOC)  GO  TO  65 
TS(J)=TJi  (1  ) 

no  Tr  i=i,Nnxi 

TJ(I)*TJ1(I) 

TJl(I)sTO 
70  CONTINUE 

41  CONf INU* 

PRINT*,"  ID  5ECt  TEMPERATURE  OISTRI9UTI0N  INTO  EARTH" 

00  7*  1=1 , NOX 1 

PPW*,"9F»TM  NOO'I  ",I,"  TFHP£RATJRE  « “,TJ(I>,“  OEG-K" 
75  CONTINUE 
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c— T<“-nr?ATUP*  C ALCULAT  I ONS  f*01  1C  TO  130  SEC- 

00  1 = 2, NOX 

O'-—.  1 CD  = 2*ALpMA*0ELTAT2/(0n."r4X(I)  *OELTAX ( I *1 ) ♦ 

CO?  l.r  6 X ( I ♦ 1)  **2) 

IF  C>rTAl  (I)  «GE«  . 5)  00  TO  202 

Oi!T12(I)=2*ALpHA*OiLTAT2/((r'"LTAX(I)  **2>  *DEITAX(I) 
r*nr_T4X <I*l>> 

TpC,rTA2(I)  ,3“,,5)  SO  TO  203 
80  COMrTMUr 

MTF.Xt  = MTpLX«-l 

1 HTEU  74/74  0°T  = 0 T RACE  FTN  4.5*41 

N1  ='|TrLX*NOTO 
00  35  J=NTFLX1,H1 
no  0 0 I = ? ,NPX 

TJK  T)=9ETAl(I)*TJ{I*l)*(l-P"T41<I)-9ETA2(I))*TJ(I) 

C*Orr A?(I) *TJ(I-1) 

(4r*“(TJl(I)-T  J1  (I*l)).LT..l)  ROTO  100 
90  COM* IS J” 

TJKM0X1)  =T  J1  (MOX ) 

lf  3 TJi( 1)  = <3*C1*(TJ<1)**4)«-C2*<3.9  3E-3  6)*<TJ<1)**1.3  33) 

C+’Ji  (?)  ) /<«»*C1*<TJ<1)  **3)  *C2*  (3.97E-05)*  <TJ<1)**.  333)  *1) 
11^  TIJr<(  = TJl(l) 

TJ1(1)=(3*C1*<TNEM  * *<»)  + C?*  (0.  9 3E-3  6)  * < TNEW  **1.333) 
0+TJ1 (2) ) /<4*C1*(TMEW  **3) *C2* (3, 97E-05) * (TNEW  **. 333) *1) 
IF  ( AOSCTJl (l)-TNEW) ,GE.TOff)  GO  TO  110 
TS(J)=TJi (1) 

00  123  1*1, NO XI 
TJ<r  ) = TJ1  (T ) 

TJK  I)  = TO 
123  OOMTisjr 
»5  COM*  TUNE 

PPTST*,”  130  SECt  TEMPERATURE  DISTRIBUTION  INTO  EARTH" 

DO  12*  1 = 1,  NO XI 

ORTNT*,"0E°TM  NODE  I “,I,“  TFSPERATJRE  = ",TJ(I),"  OEG-K” 
125  CO»riNUE 
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C TXMsr^ruor  CALCULATIONS  e-pON  190  TO  1C00  SECONDS 

00  iTC  T = ?,NOX 

or-’"  1 (I)  =2*ALPHA*0ELTAT37(OFLTAX  (I)  * DEL  TAX(l4-l)  + 
C0ELrcx<T*i)**2) 

IP(TrT(U(T).r,r.,5)  50  T0 

OETl?<I>  =2*ALP4A*0ELTAT3/( (nELTAXtl)  **2)  fOELTAX(I) 

C*Dr. TA  X ( I +1 ) ) 

IT  O ~T  A?  (I)  . SE. . 5 ) SO  TO  23? 

13S  OOMriNUE 
N?  = Nin 
N3='11«-NTTC 

no  it.;  j=n?,n? 

00  1?*  I = ?,NOX 

TJ1(  I)  = 3ctaki)*t  J(I  + l)f  (1-PETA1  <I)-nETA?(I>)*TJ(I) 

C + Srr  A-><n*TJ(I-l) 

Tc(AP>(TJl(X)-TJl(Ifl))#LT««l)  SO  TO  155 
159  COfir  TMJr 

TJ1CM0X1) =TJl (MOX ) 

155  TJ1(  i) = ( 3*C1*  (T  J(i)**4)+C2*  H.  03E-3  5)  * < T J ( 1 ) **1 , 3 33) 

C+TJ1  (2)  )/<L*Cl*  (T  JC1)  **3)fC2*  <3 . 97E- 35) *< TJ (1 )* *,  333) +1) 
16!)  TN^t'JKl) 

TJi<l)=(3*Ci* (TNEW  **4)+C2* (1. 9 3E-3 6 ) * < TNEW  **1.333) 
O+TJl  (?) ) / (4*Cl* (TNEW  **3>+02*<3.07E-05> *<TNEW  **.333>+l) 
IF(!\ns<T Jl(l)-TNEW).GC.TOF)  SO  TO  160 
TSCJ)=TJ1  (1) 

00  170  I = 1»  N0X1 
TJ(T)=TJ1(T) 

TJ1( I)=TO 
179  CONMNJE 
140  CON'TNJE 

°PT‘!T* , '*  10?C  SEC » TE^P-RATIRF  DISTRIBUTION  INTO 
C 

00  170  I*i,NOXl 

P»T3T*f-OEPTH  NODE  I ".I,**  TFNP'RATJRE  » “*TJ(I)» 

C"  n“C- <** 

175  CONTINUE 
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c frvoraftTiJPi  CALCULATIONS  F=>rK  1,909  TO  50,00  SECONOS 

oo  to*  1 = 2,  rnx 

0rTAl(I>  =?•  a Lf’H^OFLT  4T4/(0rLTA  X(I)  *nri_TAX(  I*  II 
0 + 0-LTttX<I-*-l>**?) 

IF(3FTA1  (I),r,;„5)  00  TO  295 

9rTA’(I>  =2*  ALPHA’OELTA^,/  ((OiLTAX(I)  **2>  f CELT AX (I) 
0*nr_TAX(I«-l>> 

IF  P5TA2  (T>  .0E..5)  GO  TO  207 
180  CONTINUE 
N4=M3*1 
M5=>l3f  noto 
00  t 9?  J=N4,N5 
00  ?fO  T = 2 , M 0 X 

TJK  nr^TJKT)  *T  J(IH)  M1-"ETA1  (I)  -0ETA2CI))  *TJ( I) 

Ovnrr  A?(i>  *r  j(i-i> 

Tr(*o«f(TJl(T)-TJl(I*l>),LT..l>  GOTO  205 
200  90NrTNJ- 

T Jl<  M?X  1 ) =T J1  < MTX  > 

205  TJ1C 1)  = ( 3*2 1* (T  J <1)**4) +C2* (9.D3E-36)*  <TJ(1> **1.3 33) 

C+TJ1 (?) )/ (U*Cl*<T J(l) **3) +32* (3.97£-05)*(TJ(l)**. 333) +1) 
210  TNr4  = TJKl) 

TJ1(  1)  = <3*:i*  (TNrM  **4)+C2*  n.93E-3f>)*(TNEW  **1.333) 

C*T  J1  (2)  ) / (4*0l*(TNE,(  **3)  *02*  (3.  97E-05)  * (TNEW  **.  333)+l) 
IF  (An5(Tji(i)-TMEW).G,:*T0c)  GO  TO  210 
TS(J)=TJ1  (1) 

IF  (T5( J) .LT.TMIN)  GO  TO  225 
00  ? 1 F T=i,MOXl 
TJ(T)  = TJ1(I) 

TJ1(  T)=T0 
215  CCUrIN'JE 
193  CONTINUE 

P° I*IT*, ” 53 CO 0 SFOt  TEMPERA TJ?£  9ISTRI9UT70N  INTO 
0 rft*T9“ 

00  ? 23  1*1, MO XI 

PPMr*, "DEPTH  N90E»  "fIf"  TfNP-RATJPE  = ",TJ(I),M  OEG-K" 
220  CONT TNUr 
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C-- 

22S 


fi^t  ? 5 fq  PRINT OUT  IN  .1  SEC  I NT"R VAL S 

05 JMT* 

N*>  = ,»  /DELTAT1 
N?  = . l/D^LT^l 
90  J = 1,N^>,N7 

rT=( j-i>  * n“lt a t t 

point*, et,“  srr*  t = " , ts  ( j) 

233  COM'INUE 

C 2 SEC  T9  10  SEC  PRINTOUT  IN  1 SIC  INTERVALS 

N«  =N5*2 
N9=.»7*li 

90  ’43  J = N9,  NTFLX  »M9 

ET=( J-l) ♦0ELT6T1 

POINT*, -T,"  SIC*  T = ",TS(J> 

249  CONTINUE 

C ir  to  1C9  SEC  PRINTOUT  IN  19  SEC  INTERVALS 

Nl=(lC/9ELTATl)-9 
N2=(  13/9ELT  AT  1)  +*42 
M?=( lr/9ELT  AT 1) *342 
N1  'rNT'tTU? 

00  2ED  J = N10,N1,1J 
ET=(J-  N1)*0ELTAT? 

»RINT*,rT,"  SEC*  T = ",  IS  ( J) 

259  CON^IN’Jr 

C 1C'  TEC  TO  1(1  C’ 9 SEC  °oiNTOUT  IN  190  SEC  INTE’sVALS 

N1 1=  Nl *5  9 

90  ? f> 9 J=N11,N3,S? 

ET=(J-  H2)*0CL*”ATX 

POINT*, ET,"  SEC*  T r ",TS(J> 

263  CO‘riNJE 

C l,C*f  SEC  TO  13,000  SEC  PRINTOUT  IN  1,000  SEC  INITIALS 

N12=N3*?3G 

N13=NS*ia?* 

00  27:  J = N12,N13,  2C3 

ET  = ( J-  N 3 ) * OELT AT  4 

“PINT*,  ET,"  SEC*  T = ",  TS(  J) 

270  CONTINUE 

C 1C,” r 9 SEC  to  69,903  SEC  PRTN  T?  JT  IN  5,030  SEC  INTERVALS 

N14rNlT*2C?C 
Nls=N13f 009? 

00  ?f>:  J®  Ml  4,  NIC,  1*90 

ET  = (J-  NT) * CELT  AT  4 

POINT*, ET,"  SFC  * T * ",  TS  ( J) 

269  CON'IN'JE 
r,n  ro  3?o 
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290 

291 

292 

293 

294 

, 299 

296 

297 
30  C 
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PDI«JT*, -grTfi!  c *• , I , *•  > 19  5 c * .9  FOR  DZLT4T1" 

GO  t o vo 

0<7I'JT*  , "9CT A?  CM,")  IS  GE.  .5  FOR  DELT4T1" 
GO  rO  300 

PRIMT*,"9ETA1  (H,I,-)  IS  GF • .5  FOR  D‘LT4T2” 
GO  r0  3"3 

p°I'JT*,"3FTA2  CM,”)  IS  GF.  .5  FOR  DFLTAT2” 
GO  r0  3*0 

°RT'< T* , “OFT A1  (",I,")  IS  GF«  .5  FOR  DFLTAT3" 
GO  r 0 35  0 

DoivjT*.“9FTft2  CM,”)  IS  G£.  .5  FOR  OFLT4T3“ 
GO  T0  3?9 

PPI'r*,”9ETAl  (",I,")  IS  GF  • .5  FOR  0CLT4T4" 
GO  T0  3P9 

P°  T'JT*  , ” GET  A2  (",I,")  IS  GF  • .3  FOR  0FLT4T4" 

STO3 

ENO 


Appendix  F 

Computer  Program  HTE 


HTE5  is  a short  computer  code  designed  to  calculate  the  Increase 
in  the  earth's  temperature  caused  by  a deposition  of  fission  products. 
In  the  accompanying  listing,  the  term  1.23/((.l*DELTAT*j)**1.2)  repre- 
sents the  thermal  flux  input  from  the  absorption  of  betas  by  the  earth. 
HTE5  is  similar  to  HTEA  and  uses  the  differenced  heat  transfer  equation 
for  a constant  space-time  mesh.  This  code  differs  from  HTE4  in  that 
radiation  by  the  surface  is  not  considered.  At  the  relatively  low 
temperatures  encountered,  this  omission  causes  a negligible  effect. 


Glossary!  HTE5 

ALPHA  - THERMAL  OIFFUSIVITY  <CM**2/SEC) 

BETA  DIMENSIONLESS  PARAMETER  IN  THE  HEAT 
CONDUCTION  EQUATION. 

OELTAT  - TINE  INTERVAL. 

OELTAX  - SPACE  INTERVAL. 

K - THERMAL  CONDUCTIVITY  (CAL/CM  SEC-DEG  K> 

NOT  - NUMBER  OF  TIME  INTERVALS. 

NDX  - NUMBER  OF  DEPTH  INTERVALS 

TOF  - DIFFERENCE  BETWEEN  SUCCESSIVE  ITERATIONS  OF  THE 
SURFACE  TEMPERATURE.  BELOW  WHICH  ITERATIONS 
ARE  TERMINATED 

TJ(I)  - AN  AR»AY  FOR  THF  TEMPORARY  STORAGE  OF 
THE  VALUFS  OF  T(I,JI  AT  TIME  LINE  J 

TJKII  - AN  ARRAY  FOR  THE  TEMPORARY  STORAGE  OF 

THE  VALUES  OF  T(I.Jfl)  AT  TIME  LEVFL  J»1 

TO  • SURFACE  TEMPERATURE  AT  TIME  ZERO  < EQUAL  TO  0 DES-XI 

TSCJJ  - THE  SURFACE  TEMPERATJRE  AT  TIME  LEVEL  J 
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Program  Listing t HTE< 


1T":,!'TM  T^(  n <■>  l)  , T J(?l)  , Til  { ?1) 
ao»  , al  =>ma  ,* 

-HM.  t ,,nx  ,\iOT  , - 1? 

_lt a y,r,ru tit 

q.--r*  r (AI.',M4*nrLTaT) /n“LT«V*“? 

I«*  P -T  A.  0 00  ”0  on 

t--,i  '>’■«  -\r|_rft  7 

n=  »n  v*  •>“  (_  r?  y 
mot:.  =>nr  4- 1 
•ip  y*  =‘n/  *1 

r»o T-r* , “T i'*r  r>n=*rMr^r  = ,,,n=,i  TaT,”  s“0“ 

p.->T"T»,  -)c  t T H '■  IN3?rvP*ITC  = “tMOT 

or-r.iT.  t "TOTAL  TI^r  r S " 

03T  |*»(“>)||v|n'5  Oc  rlr"’rH  - ",  *px 

OOT'IT.,  •tj^otm  T,IC3“1I*JT  S ",rL  TAX,  “ CM” 
oo  T.i-r.  ? ••.qroTq  = CM” 

nr  x-it- ,»qrTA  = ”,Q-Ta 
os  T‘!~*  , ••  ALPW  A = ”,  ALPHA,” 

os !•'*»,"<  = •*,<,••  cal <•• 

on  T-|  t * 

00  1 ' T = 1 ,Mnyi 

Tjr  ) = :. 

TJK  T>  = - . 

rnilT  Tk|'|T 

00  r j=  ? , *jr»T  l 

no  * <■  Tr?,\joy 

TJV  T)  = * m*TJ(I*T>  ♦ (1-?«P“TA)  *TJ(I)  4.0^TA*TJ(I-1) 
Tr  ('  (T)-TJi  (Tfi)  ) , LT  , , l ) 00  T0  60 

00‘|T  t.jjt 

TJK  M-tyt)  =Tjt  (vpy) 

TJ1  f M r n-LTAX/r)  Ml  .??/((.  1*0£LTAT*J)  **1,2) 
T*r<sTji.  (i) 

TJ1<  ll  = <0-LTAX/<)  M 1.  '»•*/({,  4 * AT  LTAT*  J)*  *1,2) 

O-tr  m x 5)  *T»jr,'l  *M.33>  ♦TJK'!) 


1 


* 


Tr  n °S(T  Jim-'1' l^'O  . G".T*r>  -0  TO  ss 

T<j(l)=TJl  (1  ) 

on  7'  T=l|M«yi 
tj  r )~r jt  m 
T J 1 < T)  = * . 

7 (*  T*1'J“ 

4-  ^0""  TM'I~ 

no  ■»«!  ’•» 

r>r  T‘'T* 

°3T*!T*,  •,TT‘<r  CT*<C“  O’Et>0SlTTC'r 
00  3 r J=?, 
r'T=  l*'*"LT’'T 

r>~ srn:  t = oso-tc* 

sr  coM'TMjr 

no  * i j = 4''f??C,»o 
~T  = l*0£L  Ta* 

y.A  M T.r5  74/74  O^Trl  T’ftOf 

p^T.'-e  t rr  ,*•  <?r^.  r - OET,-  If" 

SI  0PM7T’U“ 

no  s ? jr4  ;r  ,*>'  n i, 

rT-j.n'Ei.Ta'' 

PRIv't*  SPO*  T = "jTSU),"  OCG-K** 

s*>  oo*irTMJ* 

on  3 7 

ft=  j *o“l  t at 

PPT..T*  Sic*  T = D-r.-K*’ 

S3  OO'rtM'jr 

on  T • T* 

OP  l IT* 

OS  OO’OiT'O 

00  porvr*  ,"S-'rH  T<?  r,rt  #s» 

i". 1 nTn" 

"Nn 


FTN 
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Appendix  G 

DA/S  Sortie  - Neutron  Activation  Geometry 


burst  point 


Fig.  9.  DA/S  Sortie  - Neutron  Activation  Geometry 

The  following  relations  apply  to  Fig.  9. 

R2  » X{;+  HOB2 

X2  - X2+  D2-  2X  D cos(2rr-0) 
t>  a a 

B2  - X2+  D2-  2XftD  cos(2ir-G)  + HOB2 

* X2  - R2-  ALT2 

a a 

R2  - R2-  ALT2+  D2+  HOB2-  2(R^-ALT2)1^2D  cos(2tt-0) 

I 
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